Static charging of graphene and graphite slabs 
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The effect of external static charging of graphene and its flakes are investigated by using first- 
principles calculations. While the Fermi level of negatively charged graphene rises and then is quickly 
pinned by the parabolic, nearly free electron like bands, it moves down readily by removal of electrons 
from graphene. Excess charges accumulate mainly at both surfaces of graphite slab. Even more 
remarkable is that Coulomb repulsion exfoliates the graphene layers from both surfaces of positively 
charged graphite slab. The energy level structure, binding energy and and spin-polarization of 
specific adatoms adsorbed to a graphene flake can be monitored by charging. 



Graphene 1 is a semimetal having conduction and 
valance bands which cross linearly at the Fermi level 
(Ep). The resulting electron-hole symmetry reveals it- 
self in an ambipolar electric field effect, whereby under 
bias voltage the charge carriers can be tuned continuously 
between electrons and holes in significant concentrations. 
Excess electrons and holes can be also achieved through 
doping with foreign atoms. [2HI] For example, adsorbed 
alkali atoms tend to donate their valence electrons to ir*- 
bands of graphene. The excess electrons results in the 
metalization of graphene. [5] Hole doping is achieved by 
the adsorption of bismuth or antimony. [5] However, the 
system remains electrically neutral through either way of 
doping. Recently, carrier concentration and spatial dis- 
tribution of charge are also changed for very short time 
intervals by photoexcitation of electrons from the filled 
states leading to the photoexfoliation of graphite. [THS] 

In this letter, we demonstrate that the properties of 
graphene can be modified either by direct electron injec- 
tion into it or electron removal from it; namely by charg- 
ing the system externally. Remarkably, the Coulomb re- 
pulsion exfoliates the graphene layers from both surfaces 
of positively charged graphene slab. This result may be 
exploited to develop a method for intact exfoliation of 
graphene. In addition to exfoliation, the energy level 
structure, density of states, [TD] binding energies and des- 
orption of specific adatoms can be monitored by charging. 

Our results are predicted through first-principles 
plane wave calculations carried out within density func- 
tional theory (DFT) using projector-augmented wave 
potentials. [TT] The exchange correlation potential is ap- 
proximated by local-density approximation (LDA). We 
also performed GGA+vdW (generalized gradient approx- 
imation including van der Waals corrections [12J) for a 
better account of VdW interlayer interactions between 
graphite slabs. A plane-wave basis set with kinetic en- 
ergy cutoff of 500 eV is used. All atomic positions and 
lattice constants are optimized by using the conjugate 
gradient method, where the total energy and atomic 
forces arc minimized. The convergence for energy is cho- 
sen as 10 -5 eV between two steps, and the maximum 
force allowed on each atom is less than 0.01 eV/A. The 
Brillouin zone (BZ) is sampled by (15x15x5) special k- 
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Figure 1: (Color Online) Energy band structure of charged 
and neutral graphenes. (a) Positively charged graphene 
by p=+0.20 e/atom. (b) Neutral, (c) Negatively charged 
graphene by /9=-0.05 e/atom, where excess electrons start to 
occupy the surface states. Zero of energy is set to Fermi level, 
(d) Planarly averaged charge density (A) of states, \?5_io, of 
neutral graphene. (e) Charge contour plots of the lowest sur- 
face state, *l/5 = \& s in a plane perpendicular to graphene. (f) 
Same as (d) after charging with p=-0.05 e/atom. (g) Varia- 
tion of lattice constant a of graphene as a function of charging. 



points for primitive unit cell. Calculations for neutral, as 
well as charged systems are carried out by using VASP 
package. [13] Two-dimensional graphene is treated within 
periodic boundary conditions using the supercell method 
having more than 50 A separation between adjacent lay- 
ers. The amount of charging, p, is specified as either 
positive charging, i.e. electron depletion (p > 0), or neg- 
ative charging, i.e. excess electrons, in units of ± elec- 
tron (e) per carbon atom or per unit cell. For charged 
calculations, additional neutralizing background charge 
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Figure 2: (Color Online) Exfoliation of graphene layers from 
both surfaces of a 3-layer graphite slab (in AB-stacking) 
caused by electron removal. Isosurfaces of difference charge 
density, Ap, show the electron depletion. The excess charge 
on the negatively charged slab is not sufficient for exfoliation. 
The distributions of planar averaged charge density (A) per- 
pendicular to the graphene plane are shown below both for 
positive and negative charging (calculations are performed by 
GGA+vdW). 



is applied. [14 

The work function of neutral graphene is calculated 
to be 4.77 eV. Lowest two parabolic bands \l/5 and ^ 
in Fig. [jjb) have effective masses m*=1.05 and f.02 m e 
(free electron mass) in the xy-plane parallel to the atomic 
plane of graphene. Hence they are nearly free electron 
like (NFE) in 2D , but they are bound above the graphene 
plane. As shown, in Fig. [TJd) and (f), these "surface" 
states[T5] can be expressed as ^> s ~ e* U' p H $(;?), where 
r|| and ley are in the xy-plane. Parabolic bands at 
higher energies becomes NFE in 3D. When the electrons 
are removed, the Fermi level is lowered from the Dirac 
point and positively charged graphene attains metal- 
lic behavior as in Fig. [ija). At the end, the work 
function increases. However, under negative charging, 
whereby electrons are injected to the graphene, Fermi 
level raises above the Dirac point and eventually becomes 
pinned by NFE parabolic bands as in Fig. go). These 
parabolic NFE bands start to get occupied around p=- 
0.015 e/atom (or surface excess charge density er=-0.0926 
C/m 2 ). Upon charging the bound charge of ^5 states 
are further removed from graphene as shown in Fig.JIJf). 
This situation can be interpreted as the excess electrons 
start to spill out towards vacuum. Figure |TJg) shows 
another important effect of charging where the lattice 
constants increase with positive charging. On the other 
hand, negative charging has little effect on lattice con- 
stants, since the excess electrons mostly spill out. 

The effect of charging on a graphite slab consisting 
of 3 layers of graphene is better seen in Fig. [2j When 
negatively charged, the excess electrons are mainly ac- 
cumulated on both surfaces, but with smaller amount at 



the middle layer. The effect of charging on structure is 
minute, since the bonds are intact and the excess elec- 
trons rapidly spill out towards the vacuum. However, 
the situation is dramatically different for the case of pos- 
itive charging. The charge isosurfaces in Fig. [2] shows 
that positive charge, occurs mainly on both surfaces 
(i.e. first and third graphene layers), whereas the middle 
graphene has relatively small positive charge. This is an 
expected result for a metallic system. The interlayer in- 
teraction in the neutral 3-layer slab is attractive and is 
calculated to be 17 (36) meV/atom calculated by LDA 
(GGA+vdW), which becomes even weaker upon depop- 
ulation of 7r-orbitals. GGA+vdW calculations predict 
that a threshold charge, Q=0.16 e/cell gives rise to exfo- 
liation of two outermost layers. LDA calculations yield 
relatively lower threshold charge of Q=0.14 e/cell. We 
also performed a systematic analysis of exfoliation for 
thicker slabs consisting of 5-10 layers of graphene. We 
found that the threshold charge increases with increas- 
ing slab thickness. However, our analysis based on the 
planar averaged charge densities suggest that the exfoli- 
ation of outermost layers occurs when approximately the 
same amount of positive charge is accumulated on the 
outermost layers. For example, the exfoliation of 3-layer 
and 6-layer graphene flakes take place when their outer- 
most layers have positive charge of 0.065 and 0.066 e/cell, 
respectively. On the other hand, increasing of threshold 
charge by going from 3-layer to 6-layer occurs due to 
the charge spill to the inner layers. This situation can be 
explained by a simple electrostatic model, where the out- 
ermost layers of slabs is modeled by uniformly charged 
planes, which yield repulsive interaction independent of 
their separation distance, i.e. F cx q 2 /(A-eo), where q 
is excess positive charge per unit cell with the area A. 
Nonetheless, these values of charging are quite high, and 
can be attained in small flakes locally by the tip of Scan- 
ning Tunnelling Microscope. 

Ultra-fast graphene ablation was directly observed by 
means of electron crystallography. 7J Carriers excited by 
ultra-short laser pulse transfer energy to strongly cou- 
pled optical phonons. Graphite undergoes a contrac- 
tion, which is subsequently followed by an expansion 
leading eventually to laser-driven ablation. [7] Much re- 
cently, the understanding of photoexfoliation have been 
proposed, where exposure to femtosecond laser pulses has 
led to athermal exfoliation of intact graphenes.[5] Based 
on time dependent DFT calculations (TD-DFT), it is 
proposed that the femtosecond laser pulse rapidly gen- 
erates hot electron gas at ~ 20.000 K, while graphene 
layers are vibrationally cold. The hot electrons are spill 
out, leaving behind a positively charged graphite slab. 
The charge deficiency accumulated at the top and bot- 
tom surfaces lead to athermal excitation. [5] The exfolia- 
tion in static charging described in Fig.[2]is in compliance 
with the understanding of photoexcitation revealed from 
TD-DFT calculations, since the driving force which leads 
to the separation of graphenes from graphite is related 
mainly with electrostatic effects in both process. 
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Figure 3: (Color Online) Effect of charging on graphene flake consisting of 78 carbon atoms, (a) Isosurfaces of difference 
charge density Ap of positively charged, neutral and negatively charged slabs, (b) Corresponding spin-polarized energy level 
structure. Solid and continuous levels show spin up and spin down states. Red and black levels indicate filled and empty states, 
respectively. Distribution of magnetic moments at the zigzag edges are shown by insets. Zero of energy is set to Fermi level. (c) 
Variation of binding energy and net magnetic moment of specific adatoms adsorbed in two different positions, namely A-site 
and B-site indicated in (a). 



The effects of charging becomes emphasized when the 
size of graphene flake is small. In this respect, the flake 
behave like a quantum dot and hence energy level struc- 
ture is affected strongly. The flake we consider has a rect- 
angular shape and hence it consists of armchair, as well 
as zigzag edges as shown in Fig. [3j It is therefore an- 
tiferromagnetic ground state when neutral. Isosurfaces 
of difference charge density, Ap, of the same flake for 
three different charge state are shown in Fig. |3£a). Ap 
is calculated by subtracting the total charge density of 
the neutral flake from that of charged ones. For a better 
comparison, charge density of the neutral flake is calcu- 
lated using the atomic structure of the charged ones. It 
is seen that the edge states due to zigzag edges are most 
affected from charging. In Fig. [3] (a), while charge is de- 
pleted mainly from edge states, excess electrons are ac- 
cumulated predominantly at the zigzag edges. As shown 
in Fig. [3] (b), while the antiferromagnetic state of the 
flake is unaltered, charging causes emptying and filling 
of HOMO and LUMO states, changing of level spacings 
and their energies relative to vacuum level. Additionally, 
magnetic moments of zigzag edge atoms are strongly af- 
fected depending on the sign of charging in Fig. [3] (b) . In 
particular, the binding energies and magnetic moments of 
specific adatoms depends on its position and charging of 



the flake. In Fig. [3] (c) we consider Li and Ti, which nor- 
mally adsorbed to graphene by donating charge. Gener- 
ally, the binding energies increases (decreases) with posi- 
tive (negative) charging. We also found that the effects of 
monopole and dipole corrections on the effects of charg- 
ing on the binding energies is minute. For example, the 
binding energy of Li, when 2 electrons are removed, in- 
creases from 2.756 to 2.764 upon corrections. However, 
the effect of charging becomes more pronounced when the 
adatom is placed close to the edge of positively charged 
flake since the additional charges are mostly confined at 
the edges. Similarly, the magnetic moment at the adatom 
site varies depending the adsorption site and charging 
state of the flake. 

In summary, we revealed the dramatic effects of static 
and external charging of graphene and its flake. Charging 
through electron depletion of graphite surfaces leads to 
exfoliation of graphene. We also show that the binding 
energy and local magnetic moments of specific adatoms 
can be tuned by charging. 
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